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Abstract 

We present what is to our knowledge the most complete 1-D nu- 
merical analysis of the evolution and the propagation dynamics of an 
ultrashort laser pulse in a Ti:Sapphire laser oscillator. This study con- 
firms the dispersion managed model of mode-locking, and emphasizes 
the role of the Kerr nonlinearity in generating mode-locked spectra 
with a smooth and well — behaved spectral phase. A very good agree- 
ment with preliminary experimental measurements is found. 



Pulse generation in mode-locked lasers has been a subject of intense re- 
search over more than two decades. In more recent descriptions of Kerr-lens 
mode-locked (KLM) ultrafast lasers, pulse formation is assumed to rely upon 
a soliton-like mechanism, based on the master equation approximation P] , 
where the combined and balanced action of self-phase modulation (SPM) 
and group-delay dispersion (GDD) is the basis for the generation of a short 
pulse. In this picture the shortest pulse durations are obtained for cavity 
configurations with minimum net GDD. This approach gives a good descrip- 
tion of actual systems but does not take into account the pulse propagation 
dynamics in the single optical elements that comprise the laser cavity and 
usually does not include higher order dispersion terms, which are known to 
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significantly affect sub— 10 fs pulses such as those generated with today's 
state-of-the-art laser oscillators [21 El- As an alternative to the master equa- 
tion approach, the evolution of laser pulses has also been modeled numerically 
jH |3] , even if the actual sequence of intracavity components was not fully 
taken into consideration. In fact, the actual ordering of the optical elements 
is behind a new model for the generation of ultrashort pulses later introduced 
by Chen eta/. [Hj, on its hand directly related to the nonlinear propagation 
of pulses in dispersion managed communication fibers. They identify a solid 
state mode-locked laser as another example of a system where dispersion 
managed solitons (DMS) can be observed. DMS are stable soliton-like so- 
lutions of the nonlinear Schrodinger equation which are known to occur in 
optical media with a periodical change of sign in the GDD , such as a fem- 
tosecond laser cavity. The main difference from standard soliton propagation 
is that the spectrum and temporal profile of DMS periodically broaden and 
recompress as the pulse crosses regions of opposite GDD. 

Here we present what is to our knowledge the most detailed 1-D nu- 
merical simulation of a prism-dispersion controlled linear laser cavity, which 
includes build up of the laser pulse from noise through the action of active- 
and Kerr-lens mode-locking, the measured reflectivity and phase distortion 
of every optical element in the cavity, the measured gain bandwidth of the 
Ti:sapphire crystal, and the propagation inside the active medium, by nu- 
merically solving the corresponding nonlinear Schrodinger equation in the 
presence of gain. In particular, great attention is devoted to intracavity pulse 
formation and propagation, showing how the spectrum and its phase evolve 
as the pulse crosses the crystal, is reflected off the intracavity mirrors, and 
goes through the negatively dispersive prism line. This study confirms the 
dispersion managed model of mode-locking, showing that a spectrum extend- 
ing from 700 to 950 nm with a smooth and nearly flat phase can be obtained 
using commercially available ultrafast optical elements, in agreement with 
recent experimental work[8j. 



FigureH(a) shows a schematic diagram of the Ti: Sapphire laser oscillator, 
which consists of: an active crystal of length L enclosed between two focusing 
mirrors Mi and M2, a fiat folding mirror M3, an output coupler OC and 
a silver high reflector at the cavity ends, a pair of fused-silica prisms for 
dispersion compensation, and an active amplitude modulator M{t) to form 
the initial pulse from noise. The evolution of the pulse inside the cavity is 
described by the following iterative procedure which connects the spectral 
amplitude of the field, Ak+i{uj), for the (k + l) — th passage inside the cavity. 
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with the field envelope Af^it), obtained from the previous passage: 

J —oo 

Ak,{uj) = V{Ak,{uj)) 



where Ri and (pi are the refiectivity and phase of the i — th optical element 
respectively, M{t) is an initial active modulation used to start mode-locked 
operation jl], and Ak^^-j^^tu) = V{Ak^{uj)) is the spectral amplitude at the 
crystal output, obtained by numerically solving the following propagation 
equation inside the crystal: 

dA{z,t) 
dz 



{il3{iu) + g{u))A{z,u)e~''^'duj 
+^J\A{z,t)\'A{z,t) (2) 



with ^4(0,^) = Ak^{t). Here (3{u) and g{uj) are the crystal phase distor- 
tion and gain profile per unit length, and 7 = 2.2 x 10~^W^~^ cm~^ is 
the estimated nonlinearity coefficient. Moreover, self- amplitude modulation 
(SAM) induced by the Kerr nonlinearity is included as a nonlinear inten- 
sity discriminator in the time domain^, here modeled as a supergaussian: 
K{A) = exj9[-|(^)™], with Pq = max{\A\'^), m = 24 and a = 0.47. 
These parameters were determined by trial-and-error, even though the final 
spectrum will not significantly depend on the exact values, provided that 
mode-locking operation is established. The crystal gain and the mirror re- 
fiectivities and phase distortions were directly obtained from measurements 
performed on actual optical components, which comprise a relatively thick 
(4.5mm) Ti:Sapphire crystal (Crystal Systems Inc.), commercially available 
standard ultrafast laser mirrors designed for 850 nm (Spectra-Physics and 
TecOptics) and a 3.5% output-coupler (Spectra-Physics). Figures[T](b) and 
(c) show the total GDD and the total cavity gain Gtot for one round trip de- 
fined as GtotiX) = RiiXfR2iX)^R3iXfRocWGiX)^ with G(A) = e^^^ the 
total crystal gain normalized at 1.04. As initial condition, a low amplitude 
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random noise with flat spectral phase is assumed. The evolution of the in- 
tracavity spectral profile and phase of the pulse prior to entering the output 
coupler are shown in figureEl after (a) 20, (b) 500, and (c) 1000 round-trips of 
the pulse inside the cavity. After the first 20 round-trips, the spectral phase 
of the intracavity radiation shows strong fluctuations due to the accumulated 
linear phase distortion of the optics, as expected from figuren(b). Then a 
spectrum with increasingly smooth spectral phase centered around 850 nm 
builds up from noise mainly through the action of the crystal gain, the ac- 
tive modulation and SAM, resulting in the spectrum shown in figure|21(b) 
at a pulse peak power of 3.5 x 10^ W. As the pulse intensity increases, new 
mode-locked frequencies are generated via SPM and the spectrum is further 
broadened until a steady-state is reached at a pulse peak power of 5 x 10^ W 
(see figure|21(c)), as a result of the interplay between SPM, phase distortion 
and the finite bandwidth of the total cavity gain (see figure[T](c)). To il- 
lustrate the validity of the simulation code the measured and the simulated 
spectra obtained after the OC are given in figure|21(d), revealing a very good 
agreement between experimental measurements and theoretical predictions. 
The steady-state spectrum of figure|21(c) is reproduced every round-trip but 
it has a different spectral profile in different dispersion regions of the cavity. 
FigureEl shows how the cavity can be put in analogy with a dispersion man- 
aged fiber made of a region of nonlinear propagation and positive GDD (the 
active medium crossed two times), and a region of linear propagation and 
negative GDD (the prism-pair crossed two times). The OC and the silver 
mirror are placed in the middle point of each dispersion region. In correspon- 
dence with the points A, B,C, D, E and F of figureEl figurelH shows how the 
steady-state pulse spectrum, phase, and temporal profile evolve when cross- 
ing the crystal, going through the dispersion compensating prism-pair, and 
coming back to the OC. As the pulse propagates towards the silver mirror, 
it enters the active medium with a positive chirp (point A in figures^] and|3]) 
and its spectrum broadens while its temporal width increases due to the 
concomitant action of positive GDD and SPM (points B and C). It can 
be observed that spectral broadening mainly takes place in the first half of 
the crystal (from A to B) since temporal broadening due to GDD decreases 
the strength of SPM as the pulse penetrates more into the crystal. When 
the pulse crosses the two prisms and reaches the silver mirror (point D), 
the acquired positive chirp is partially compensated for, while the spectrum 
remains unaffected, resulting in a shorter and more intense pulse. This is 
the point of maximum intracavity spectral width, minimum pulse duration 
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and maximum peak intensity. In contrast to the case of point A, when the 
pulse enters the crystal on its way back to the OC (point E), it is negatively 
chirped by a second passage in the prism sequence: nonlinear propagation 
now results in spectral narrowing and phase flattening until the pulse is trans- 
form limited at a depth of 3 mm within the crystal (point F). This is the 
point of minimum spectral width. As the pulse continues its way to the OC 
and goes back to the silver mirror, its spectrum is broadened, recovering its 
maximum width at points C and D. The pulse behaviour described here 
confirms the dispersion managed model of mode- locking :6J, by showing the 
typical spectral and temporal breathing of DMS, and also explains the ex- 
perimental fact that spectra are broader when taken in the dispersive end of 
KLM lasers 0. Furthermore figurelH shows that the pulse propagating in the 
cavity, while varying its spectrum and temporal profile, mantains a smooth 
spectral phase due to the joint action of SAM and SPM, which are able to 
wash out from the spectral phase the modulations inherent in the net intra- 
cavity GDD. This is in agreement with recent experimental results, where 
it was shown that operation of a broadband Ti: Sapphire laser under strong 
Kerr-lens mode-locking conditions resulted in a smoothing of the spectral 
phase [3 El- Once extracted the pulse from the cavity, such a well — behaved 
phase appears to be suitable for further extracavity pulse compression, giving 
for the spectrum of point D a Fourier transform limited pulse duration of 
11 fs, also in good agreement with recent experimental results [S]. 

In conclusion we proved numerically the dispersion managed model of 
mode-locking using what is to our knowledge the most complete 1-D sim- 
ulation of an actual ultrafast Ti:Sapphire laser. Moreover, we show how 
the pulse acquires and preserves a smooth spectral phase in its propaga- 
tion inside the cavity. These results are in very good agreement with ex- 
perimental measurements of the spectrum outside the OC. A detailed ex- 
perimental demonstration of pulse evolution as predicted by this model is 
presently under development. This work was partly supported by FCT Grant 
No. POCTI/FIS/48709/2002, Portuguese Ministry of Science, co-financed by 
FEDER. M. V. Tognetti's e-mail address is marco.tognetti@fc.up.pt. 
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List of Figure Captions 



Fig. 1 (a) Schematic diagram of the Ti:sapphire laser oscillator, (b) Total 
round trip intracavity GDD. (c) Total cavity gain Gtot assuming a peak 
crystal gain of 1.04. 

Fig. 2 Normalized intracavity spectrum (solid line) and phase (dotted line) 
before the OC after (a) 20, (6) 500, (c) 1000 round trips inside the cavity, (d) 
Simulated (solid line) and measured (dashed line) final spectra outside the 
OC. 

Fig. 3 Laser cavity dispersion map. A, B, C, D, E and F are reference points 
corresponding to the spectra, phases, and temporal profiles reported in fig- 
urelH (see the text for more details). 

Fig. 4 Spectra, phases and temporal profiles corresponding to the reference 
points reported in figurelHl 
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Figure 1: (a) Schematic diagram of the Ti:sapphire laser oscillator, (b) 
Total round trip intracavity GDD. (c) Total cavity gain assuming a 
peak crystal gain of 1.04. 
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Figure 2: Normalized intracavity spectrum (solid line) and phase (dotted 
line) before the OC after (a) 20, (6) 500, (c) 1000 round trips inside the cavity, 
(d) Simulated (solid line) and measured (dashed line) final spectra outside 
the OC. 
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Figure 3: Laser cavity dispersion map. A, B, C, D, E and F are reference 
points corresponding to the spectra, phases, and temporal profiles reported 
in figurelU (see the text for more details). 
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Figure 4: Spectra, phases and temporal profiles corresponding to the refer- 
ence points reported in figureEl 
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